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Abstract 

In view of the fact that the data on neutrino mixing are still compatible with a 
situation where Bimaximal mixing is valid in first approximation and it is then 
corrected by terms of 0{\c) (with Ac being the Cabibbo angle), arising from the 
diagonalization of the charged lepton masses, such that 5 sin 2 6±2 ~ sin #13 ~ O(Xc) 
while 5 sin 2 #23 ~ O(Ap), we construct a model based on the discrete group S4 where 
those properties are naturally realized. The model is supersymmetric in 4-dimensions 
and the complete flavour group is £4 x Z4 x U(1)fn, which also allows to reproduce 
the hierarchy of the charged lepton spectrum. The only fine tuning needed in the 
model is to reproduce the small observed value of r = Am 2 un /Am 2 tm . Once the 
relevant parameters are set to accommodate r then the spectrum of light neutrinos 
shows a moderate normal hierarchy and is compatible, within large ambiguities, with 
the constraints from leptogenesis as an explanation of the baryon asymmetry in the 
Universe. 
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1 Introduction 



It is an experimental fact [1-3] that within measurement errors the observed neutrino 
mixing matrix [4] is compatible with the so called Tri-Bimaximal (TB) form [5]. The best 
measured neutrino mixing angle 9 12 is just about la below the TB value tan 2 ^ = 1/2, 
while the other two angles are well inside the la interval [1]. In a series of papers [6-9] 
it has been pointed out that a broken flavour symmetry based on the discrete group A4 
appears to be particularly suitable to reproduce this specific mixing pattern as a first 
approximation. Other solutions based on alternative discrete or continuous flavour groups 
have also been considered [10-12], but the A4 models have a very economical and attractive 
structure, e.g. in terms of group representations and of field content Q In most of the 
models A± is accompanied by additional flavour symmetries, either discrete like or 
continuous like U(l), which are necessary to eliminate unwanted couplings, to ensure the 
needed vacuum alignment and to reproduce the observed mass hierarchies. Given the set 
of flavour symmetries and having specified the field content, the non leading corrections to 
the TB mixing arising from loop effects and higher dimensional operators can be evaluated 
in a well defined expansion. In the absence of specific dynamical tricks, in a generic model, 
all the three mixing angles receive corrections of the same order of magnitude. Since the 
experimentally allowed departures of 9 12 from the TB value sin 2 9\2 = 1/3 are small, at 
most of O(Xq), with Ac the Cabibbo angle, it follows that both 9\ 3 and the deviation 
of #23 from the maximal value are expected in these models to also be at most of 0(X 2 C ) 
(note that Ac is a convenient hierarchy parameter not only for quarks but also in the 
charged lepton sector with m^/m T ~ 0.06 ~ X^ and m e /m^ ~ 0.005 ~ ^c~ 4 )- A va l ue 
of #13 ~ 0(Xq) is within the sensitivity of the experiments which are now in preparation 
and will take data in the near future. However, the present data do not exclude a larger 
value for #13 with #13 ~ 0(Xc)- In fact, two recent analysis of the available data lead to 
sin 2 # 13 = 0.016 ±0.010 at la [2] and sin 2 13 = O.OlOjftou at 1(7 [ 3 L which are compatible 
with both options. If experimentally it is found that #13 is near its present upper bound, 
this could be interpreted as an indication that the agreement with the TB mixing is 
accidental. Then a scheme where instead the Bimaximal (BM) mixing is the correct 
first approximation modified by terms of 0{Xc) could be relevant. This is in line with 
the well known empirical observation that 9 12 + Ac ~ vr/4, a relation known as quark- 
lepton complementarity [14], or similarly 9 12 + ^m^/m T ~ 7r/4. No compelling model 

Recently, in ref. [13], the claim was made that, in order to obtain the TB mixing "without fine tuning" , 
the unique finite group must be S4 or a larger group containing S4. For us this claim is not well grounded 
being based on an abstract mathematical criterium for a natural model. We do not share the definition 
of ref. [13], based on group theory, of a natural model. For us a physical field theory model is natural if 
the interesting results are obtained from a lagrangian that is the most general given the stated symmetry 
and the specified representation content for the flavons. For example, in [7-9], the authors obtain from 
A4 (which is a subgroup of S4) a natural (in our sense) model for the TB mixing by simply not including 
symmetry breaking flavons transforming like the 1' and the 1" representations of A4 (a restriction not 
allowed by the rules specified in ref. [13]). Rather, for naturalness we require that additional physical 
properties like the hierarchy of charged lepton masses also follow from the assumed symmetry and are 
not obtained by fine tuning parameters: for this actually A4 can be more effective than S4 because it 
possesses three different singlet representations 1,1' and 1". 
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leading without parameter fixing to the exact complementary relation has been produced 
so far. Probably the exact complementarity relation is to be replaced with something like 
612 + O(Xc) ~ (which we could call "weak" complementarity). 

In the present paper we use the expertise recently acquired with non Abelian finite 
flavour groups to construct a model based on the permutation group S4 which naturally 
leads to the BM mixing at the leading level. We adopt a supersymmetric formulation of 
the model in 4 space-time dimensions. The complete flavour group is S4XZ4X U(1)fn- In 
the lowest approximation, the charged leptons are diagonal and hierarchical and the light 
neutrino mass matrix, after see-saw, leads to the exact BM mixing. The model is built 
in such a way that the dominant corrections to the BM mixing, from higher dimensional 
operators in the superpotential, only arise from the charged lepton sector at Next-to-the- 
Leading-Order (NLO) and naturally inherit Ac as the relevant expansion parameter. As 
a result the mixing angles deviate from the BM values by terms of O(Xc) (at most), and 
weak complementarity holds. A crucial feature of the model is that only 6* 12 and #13 are 
corrected by terms of 0(Xc) while 6*23 is unchanged at this order (which is essential to 
make the model agree with the present data). We finally discuss the phenomenology of 
this model and the outlook. 

The article is organized as follows. In the next section we recall the general properties 
of the neutrino mass matrix leading to the BM mixing in the flavour basis. In section 3 
we explain why £4 is a simple candidate group for the task. In section 4 we discuss the 
structure of the model at the leading order (LO), where the resulting lepton mixing is 
exactly the BM. In section 5 we justify the choice of the vacuum assumed in the previous 
section, by minimizing the scalar potential of the theory in the supersymmetric limit. In 
section 6 we analyze the higher-order corrections needed to modify the BM mixing in 
order to derive an acceptable lepton mixing pattern. Section 7 contains a discussion of 
the phenomenological predictions of the model also including the constraints coming from 
leptogenesis and, finally, section 8 is devoted to our conclusion. 

2 Bimaximal Mixing 

We start by recalling that, in the basis where charged leptons are diagonal, the most 
general mass matrix which corresponds to #13 = and #23 maximal is of the form [15] 



Note that this matrix is symmetric under the 2-3 or the /j, — r exchange. For #13 = there 
is no CP violation, and phases are only of Majorana type. If we disregard them, we can 
restrict our considerations to real parameters. There are four of them in eq. Q which 
correspond to three mass eigenvalues and one remaining mixing angle, 6 I 12 . In particular, 
#12 is given by: 




(1) 



sin 2 29 




(2) 



(x — w — z) 2 + 8y 2 
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In the BM case sin 2 29 12 = 1 is also fixed and an additional parameter, for example w, can 
be eliminated, leading to: 

fx y y \ 
m vB M = \ y z x - z , (3) 
\y x — z z ) 

This mass matrix is exactly diagonalized by the BM pattern, which we can write, in a 
particular phase convention, as: 
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(4) 



In the BM scheme tan 2 #12 = 1, to be compared with the latest experimental determination: 
tan 2 9 12 = 0.45 ±0.04 (at la) [1-3], so that a rather large non leading correction is needed, 
as already mentioned in the Introduction. 

Both the TB and the BM mixing matrices are of a form that suggests that mixing 
angles are independent of mass ratios (while for quarks and charged leptons relations like 
\ 2 C ~ md/m s or m^/rrir are typical). In fact, in the basis where charged lepton masses are 
diagonal, the effective neutrino mass matrix in the BM case is given by 

m uB M = U B Mdiag(m 1 ,m 2 ,m 3 )Ul M 



where 
Mi = 




M< 



■y/2 
1 
1 




(5) 










)■ 
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-:■) 




-1 





(6) 

The eigenvalues of m v are m 1; m 2 , m 3 with eigenvectors (y/2, 1, l)/2, (—y/2, 1, l)/2 and 
(0, 1, — l)/y/2, respectively. In terms of the parameters x, y and z of eq. (J3J) we have 
mi = % + y/2y, mi = x — y/2y, 7713 = 2z — x. Clearly, all type of hierarchies among neutrino 
masses can be accommodated. The smallness of the ratio r = Am 2 nn /Am 2 4m = (|m 2 | 2 — 
|mi| 2 )/||m3| 2 — |m2| 2 | requires either \xy\ <C \z 2 \ (normal hierarchy) or \x\ ~ \z\ \y\ 
(inverse hierarchy) or \y\ <C \x\ ~ \z\ (approximate degeneracy except for x ~ 2z). From 
the general expression of the eigenvectors one immediately has a confirmation that this 
mass matrix, independent of the values of rrii, leads to the BM mixing matrix. 



3 A convenient presentation of the group 04 

We first present the argument to show that £4, the permutation group of 4 elements, 
is a good candidate for a flavour symmetry to realize the BM mixing. The group S4 has 
24 transformations and 5 irreducible representations, which are 3, 3', 2, 1 and V. 
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Table 1: Character table of S4. C, are the conjugacy classes, n the number of 
elements in each class. 



The character table of the group is shown in Table [T] and, in terms of two operators S 
and T satisfying to 

T 4 = S 2 = (ST) 3 = (TS) 3 = 1 , (7) 

all the 24 S4 transformations can be obtained by taking suitable products. Divided into 
the 5 equivalence classes (as many as the inequivalent irreducible representations) the 
transformations are given in Table [2} Different presentations of the S4 group have been 
discussed in the recent literature [12,16,17] and the choice of generators that we adopt in 
this paper is related to other existing choices by unitary transformations. 





Tr 


Products 




3 


1 


c 2 
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T 2 , ST 2 S,(ST 2 ) 2 
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ST, TS, (ST) 2 , (TS) 2 , T 2 ST, TST 2 , T 3 ST 2 , T 2 ST 3 


c 4 
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S, T 3 ST, TST 3 , T 2 ST 2 , ST 2 ST, TST 2 S 


c 5 


-1 


T, T 3 , ST 2 , T 2 S, STS, TST 



Table 2: The 24 S4 transformations obtained as products of S and T and 
divided into the 5 equivalence classes. In the second column the trace, in the 
representation 3, of all the representative matrices in each class. 



Explicit forms of S and T in each of the irreducible representations can be simply 
obtained. In the representation 1 we have T = 1 and S = 1, while T = — 1 and S = — 1 
in 1'. In the representation 2 we have: 



T 



1 
-1 



S 



-1 \/3 
1 



For the representation 3, the generators are: 
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(9) 
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In the representation 3' the generators S and T are simply opposite in sign with respect 
to those in the 3 (a glance to Table 2 shows that all elements in the equivalence classes 
C4 and C5 change sign going from 3 to 3' while those in classes Ci, C 2 and C3 remain the 
same) . 

We now recall the multiplication table for S± (R stands for any representation). The 
Clebsch-Gordan coefficients in our basis are collected in the Appendix A. 

1®R=R®1=R 
1'© V = 1 
1/® 2 = 2 
V © 3 = 3' 
1' © 3' = 3 



2®2 = 1© V ©2 
2 © 3 = 3 © 3' 
2 © 3' = 3 © 3' 



3©3 = 3' ©3' = 1©2©3©3' 
3 © 3' = V © 2 © 3 © 3' 

This description of the group S4 is particularly suitable for our purposes because the 
general neutrino mass matrix corresponding to the BM mixing, in the basis where charged 
leptons are diagonal, given by eq. ([3]), is left invariant by the unitary, real and symmetric 
matrix S in eq. (pi): 

m vB M = Sm uBM S . (10) 

More precisely, the matrix m V BM can be completely characterized by the requirement of 
being invariant under the action of S, eqs. (j9])-(10), and under the action of U, also a 
unitary, real and symmetric matrix 




m vBM = Um uBM U , £7=001. (11) 



The most general solution to the eqs. (10) and (11), is given by the matrix of eq. ([3]). 
The matrices S and U commute and generate a discrete group of the type Z 2 x Z 2 . As 
we shall see, in our model the invariance under U arises accidentally, as a consequence of 
the specific field content and is limited to the contribution of the dominant terms to the 
neutrino mass matrix. For this reason we do not need to include the U generator in the 
flavour symmetry group. 

Charged leptons must be diagonal in the basis where m v has the BM form. A diagonal 
matrix ml mi with generic entries is invariant under T given in eq. [9J 

ml mi = T + mlm{r (12) 

and, conversely, the most general hermitian matrix invariant under T is diagonal [^} By 
starting from a flavour symmetry group containing 5*4, we will realize a special vacuum 

2 This property remains true also when T is replaced by rjT where rj represents an arbitrary phase. 
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alignment such that, at LO, the residual symmetry in the neutrino sector will be that 
generated by S and U, while in the charged lepton sector will be, up to a phase, that 
generated by T. Then, by construction, the LO lepton mixing in this model will be of the 
BM type. Furthermore, a realistic model will be obtained by adding suitable subleading 
corrections to this zeroth order approximation. 

4 The structure of the model 

We discuss here the general properties of our model, which leads to the BM mixing in 
first approximation and is formulated in terms of the S4 realization introduced above. We 
formulate our model in the framework of the see-saw mechanism (even though it would 
also be possible to build a version where light neutrino masses are directly described by 
a single set of higher dimensional operators, violating the total lepton number by two 
units). For this we assign the 3 generations of left-handed (LH) lepton doublets / and 
of right-handed (RH) neutrinos v c to two triplets 3, while the RH charged leptons e c , 
[f and r c transform as 1, 1' and 1, respectively. The S4 symmetry is then broken by 
suitable triplet flavons. All the flavon fields are singlets under the Standard Model gauge 
group. Additional symmetries are needed, in general, to prevent unwanted couplings and 
to obtain a natural hierarchy among m e , and m T . In our model, the complete flavour 
symmetry is S4 x Z4 x U(1)fn- A flavon 9, carrying a negative unit of the U(1)fn charge 
F, acquires a vacuum expectation value (VEV) and breaks U(1)fn- 111 view of a possible 
GUT extension of the model at a later stage, we adopt a supersymmetric context, so that 
two Higgs doublets h Ut d, invariant under S4, are present in the model. A U(1)r symmetry 
related to R-parity and the presence of driving fields in the flavon superpotential are 
common features of supersymmetric formulations. Supersymmetry also helps producing 
and maintaining the hierarchy (h u ^) — ^«,d A where A is the cut-off scale of the theory. 

The fields in the model and their classification under the symmetry are summarized in 
Table |3j The complete superpotential can be written as 

w = wi + w v + w d , (13) 

with (we indicate with (. . .) the singlet 1, with (. . .)' the singlet 1' and with (. . 
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Table 3: Transformation properties of all the fields. 
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(R = 2, 3, 3') the representation R) 

y W e 2 y ( e ] 9 2 yf ] 9 2 

+ + (14) 

w„ = y{v c l) + MA(z/V c ) + a{v c v c iv) + b(v c v c <p v ) + . . . (15) 

w d = M v A(<pl<p v ) + gx (<pl(<p„<Pv) 3 ) + 92 (<fl<Pv) + 

+ Mf A 2 £° + + gsCUu + g&iWv) + 

+ /i (^°(^,) 2 ) + / 2 W(xiXi) 2 ) + /a W(viXi) 2 ) + 

+ /4W(Mi)3')+.... (16) 

To keep our formulae compact, we use a two-component notation for the fermion fields 
and omit to write the Higgs fields h u ^- For instance y T r c (l(pi)/ A stands for y r r c {lipi)h d / A, 
y(v°l) stands for y(u c l)h u . The powers of the cutoff A also take into account the presence 
of the omitted Higgs fields. Note that the parameters M, M^, and defined above 
are dimensionless. In the above expression for the superpotential w, only the lowest order 
operators in an expansion in powers of 1/A are explicitly shown. Dots stand for higher 
dimensional operators that will be discussed later on. The stated symmetries ensure that, 
for the leading terms, the flavons that appear in wi cannot contribute to w u and viceversa. 

We will show in Sect. 5 that the potential corresponding to Wd possesses an isolated 
minimum for the following VEV configuration: 



(fi) _ [ 1 1 A (Xi) 





(17) 



= D (18) 



A V -i J A 

where the factors A, B, C, D should obey to the relations: 

V3frA 2 + V3f 2 B 2 + f 3 AB = (19) 

D = _M* c 2 = 9lM ' + 93M * ' 92M ^ (20) 

92 2g 2 g 4 

Notice the existence of a flat direction related to an arbitrary, common rescaling of A 
and B: if we indicate with m 2 pi and m 2 Xi the soft masses of the two flavons tpi and xu we 
can assume m^,m^ < and then (ipi) and (xi) slide to a large scale, which we assume 
to be eventually stabilized by one-loop radiative corrections, fixing in this way A and 
B. In the discrete component S 4 x Z 4 of the full flavour group we can choose generators 
(S,T,i), where the imaginary unit % denotes the generator of the Z 4 factor. The flavons 
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£„ and <p v are invariant under Z$ and their VEV's are eigenvectors of the generator S 
corresponding to the eigenvalue 1, so that the corresponding breaking of S4 x Z4 preserves 
the subgroup G u generated by (S, i). In the charged lepton sector S4 x Z 4 is broken down 
to the subgroup G[, generated by the product %T . Indeed the generator iT is given by 
±diag(— i, 1, — 1), with the plus (minus) sign for the 3 (3') S4 representation. Such a 



generator, acting in the appropriate representation on the VEV's of eq. (17), leaves them 
invariant. It is precisely the mismatch, present at the LO, between the subgroups G v and 
Gi preserved in the neutrino and charged lepton sectors, respectively, that produces the 
BM lepton mixing, as we will explicitly see in this section. 

Similarly, the Froggatt-Nielsen flavon 9 gets a VEV, determined by the D-term asso- 
ciated to the local U{1)fn symmetry (see sect 5), and it is denoted by 

f=«. PI) 
With this VEV's configuration, the charged lepton mass matrix is diagonal 

v d (22) 



mi 



( (y ( ( } ) B 2 -y { J 2) A 2 + y ( J i) AB)t 2 \ 

y„Bt 

V y T A J 



so that at LO there is no contribution to the lepton mixing matrix from the diagonalization 
of charged lepton masses. In the neutrino sector for the Dirac and RH Majorana matrices 
we have 

/ 1 \ / 2M + 2aD -2bC -2bC \ 

m » = ( 1 \yv u M N = I -2bC 2M + 2aD A . 

\ 1 / \ -2bC 2M + 2aD / 

(23) 

The matrix can be diagonalized by the BM mixing matrix Ubm, which represents the 
full lepton mixing at the LO, and the eigenvalues are 

Mi = 2\M + aD- V2bC\A M 2 = 2\M + aD + V2bC\A M 3 = 2\M + aD\A . (24) 

After see-saw, since the Dirac neutrino mass matrix commutes with and its square 
is a matrix proportional to unity, the light neutrino Majorana mass matrix, given by the 
see-saw relation m u = —(m®) T M^m® is also diagonalized by the BM mixing matrix and 
the eigenvalues are 

1 1 \y 2 K 1 1 1 \y 2 \vl 1 . . \y 2 \vl 1 

\mi\ — 7= r W2 = 7= r \n^v = —, ;~r ■ 

2\M + aD - V2bC\ A 2|M + aD + \[2bC\ A 2|M + aD\ A 

(25) 

The light neutrino mass matrix depends on only 2 effective parameters, at LO, indeed the 
terms M and aD enter the mass matrix in the combination F = M + aD. The coefficients 
ye\ y^i 2/t> 2/) a an d b are all expected to be of 0(1). A priori M could be of C(l), 
corresponding to a RH neutrino Majorana mass of O(A), but, actually, we will see that it 
must be of the same order as C and D. 
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We expect a common order of magnitude for the VEV's (scaled by the cutoff A): 

A ~ B ~ v , C ~ D ~ «' . (26) 

However, due to the different minimization conditions that determine (A, B) and (C,D), 
we may tolerate a moderate hierarchy between v and i>'. Similarly the order of magnitude 
of t is in principle unrelated to those of v and v' (see section 5). It is possible to estimate 
the values of v and t by looking at the mass ratios of charged leptons: 



0.06 — ^ ~ t 



T / exp T 

— ~ 0.005 — ~ ut (27) 

"We*? "V 



m P \ rrif, 



0.0003 — ~ vt 2 . 



In order to fit these relations, approximately we must have t ~ 0.06 and v ~ 0.08 (modulo 
coefficients of 0(1)). 



To summarize, at LO we have diagonal and hierarchical charged leptons together with 
the exact BM mixing for neutrinos. It is clear that substantial NLO corrections are needed 
to bring the model to agree with the data on 6\2- A crucial feature of our model is that 
the neutrino sector flavons ip u and £ u are invariant under Z 4 which is not the case for the 
charged lepton sector flavons tpi and xi- The consequence is that tp u and £ u can contribute 
at NLO to the corrections in the charged lepton sector, while at NLO (pi and xi cannot 
modify the neutrino sector couplings. As a results the dominant genuine corrections to 
the BM mixing only occur at NLO through the diagonalization of the charged leptons. 
We will discuss the NLO corrections in section 6 after having proven that the needed VEV 
alignment is in fact realized at LO. 



4.1 The light neutrino spectrum and the value of r 

We now discuss the constraints on the parameters of the neutrino mass matrix in order 
to get the correct value for the ratio r = Am^/Am^ = (|m 2 | 2 — |^i| 2 )/||?ri3| 2 — |^2| 2 |, 
which is r = 0.032^q;qq5 at 3cr's. Like for models based on A 4 , also in this case some fine 
tuning is needed to accommodate the value of r. In fact the triplet assignment for LH 
lepton doublets and for RH neutrinos tends to lead to r ~ 1. We find useful to begin the 
presentation by analyzing the LO terms, even though a more complete phenomenological 
discussion with the inclusion of the NLO contributions will be illustrated in section 7. 



We redefine the parameters in eqs. (25) as follows: 

F = M + aD, Y = -V2bC, (28) 

so that 

Y ~ v' . (29) 
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while F, like M, a priori could be larger, of 0(1). We make the phases of F and Y explicit 
by setting 

F -> Fe i4,F Y -> IV 05 " , (30) 

where now F and y are real and positive parameters. Defining the phase difference 
A = (fry — <Pf we can explicitly write the absolute values of the neutrino masses as 

1 \y 2 \vl 



mi 



\m 2 \ 



y/F 2 + Y 2 + 2FY cos A 2A 
y/F 2 + Y 2 — 2FY cos A 2A 



1 Iw 2 ^ 2 

Kl - (3D 

Note that the phase A is related by a non-trivial relation to the Majorana CP phase 2a 2 i, 
which, by definition, is the phase difference between the complex eigenvalues ra\ and m 2 
of the neutrino mass matrix. Furthermore cos A must be positive in order to guarantee 
\ni2\ > |toi|. By defining F/Y = f, we can write the expression of the ratio r: 

4/ 3 cos A 



(/ 2 + l + 2/cosA)(l-2/cosA) ' 



(32) 



In order to have r small either we take / small or cos A small (or both). If F ~ 0(1) 
then / ~ 0(l/v'), r ~ 4/ cos A and cos A must be extremely small : cos A ~ v'r/A ~ 10~ 3 . 
We prefer to take / small, such that 

r~4/ 3 cosA. (33) 

If so, in order to accommodate the value of r, we only need, for example, 4 cos A ~ 1 and 
/ ~ 1/3. In conclusion, we have to take F ~ M ~ 0(v') and / = F/Y moderately small. 
We interpret the relation M ~ F ~ v', needed to reproduce the value of r, as related to 
the fact that the RH neutrino Majorana M mass must empirically be smaller than the 
cut-off A of the theory. In the context of a grand unified theory this corresponds to the 
requirement that M is of 0(M GUT ) rather than of 0(M P i anck ). 

With these positions, the neutrino spectrum shows a moderate normal hierarchy, with 

Ka| ~ \ - O (1) (34) 

Kl ~ o(^-] , (35) 



jy 

in units of |y 2 |f 2 /2A. At the leading order an inverted ordering of the neutrino masses is 



forbidden, as we can see from eq. (31), which, for \m2\ > |?tii|, always implies \m^\ > \m\\. 



At the LO the lightest neutrino mass |mi| has a lower bound of order 0.01 eV. Indeed, 
the only possible way to decrease \m\\ is to take Y as large as possible. By expanding eqs. 



(31) in powers of /, we have 

Kl « f\m 3 \ > fV&rnltm > 0-01 eV . (36) 
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As we will see in section 7, this lower bound can be evaded by including NLO corrections, 
but values of \m\\ much smaller than 0.01 eV would require an additional tuning of the 
parameters. 

4.2 Effective Terms 

Since the neutrino sector is not charged under the Z 4 symmetry, we have operators of 
dimension 5 which contribute to the neutrino masses and may corresponds to some heavy 
exchange other than the RH neutrinos v c . These effective terms are 

M' a' b' 

w e J f = —(lh u lh u ) + j^(lhjh u )£ u + j-^(lh u lh u <p„) . (37) 

The choice of the names M', a' and b' is not accidental, in fact the resulting mass matrix 



has a flavour structure identical to that one of M N in eq.(23): 



2M' + 2a'D -2b'C -2b'C , , 

nfjf = I -2b'C 2M' + 2a! D ^ . (3* 



-2b' C 2M' + 2a' D 



A 



As a consequence, also this mass matrix can be diagonalized by the BM pattern. The 
complete neutrino mass matrix is then given by the sum of two contributions, 

m u = -[m D v ) T M- N l m D v + , (39) 

and can still be diagonalized by the the BM scheme. Moreover, since both these terms 
have the same structure, m v can be parameterized by only three parameters: in principle 
we can absorb the effective contributions in the see-saw ones. 
If we consider y ~ 0(1), F ~ v' and M' ~ 0(1) we get that 

m li2 = 0(i|) + 0(|) (40) 
^-0(± V f) + 0( V f). (41) 

If also M' ~ 0(v') then the effective terms are even more suppressed. In conclusion the 
effective terms are subdominant for the interesting domain of parameters. 



5 Vacuum alignment 



In this section we show that the superpotential Wd given in eq. (16) has an isolated 



minimum that corresponds to the VEV's in eqs. (17) and (18). At the LO, the equations 



for the minimum of the potential can be divided into two decoupled parts: one for the 
neutrino sector and one for the charged lepton sector. Given this separation, for shorthand 
we can omit the indexes I and v for flavons and driving fields. The driving fields are 
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assumed to have vanishing VEV's. In the neutrino sector, the equations for the vanishing 
of the derivatives of Wd with respect to each component of the driving fields are: 



ilWx + 9l {ifl -<£) + g 2 &i =0 (42) 
M^A<p 3 - 2g m <p 2 + g 2 & 3 = (43) 
M V K V2 + 2g m <p 3 + g 2 £<p 2 = (44) 



M 2 A 2 + + gsC + gtkpl + 2<p 2 <p 3 ) = . 



(45) 



A solution to these equations is given in eqs. (18), (20). The triplet VEV is eigenvector 



of the operator S, which in the 3 representation is given by the matrix in eq. Q, with 
eigenvalue 1. 

For the charged lepton sector the equations are: 



h(fi - ws) + Mx! - X2X3) + ~y h(^2X2 + ^sxs) 







-<PiXi + 2 (^2X3 + ^3X2) 



U(V3X3 - ¥2X2) = 

U(.-<PiXa ~ V2X1) = 
h(<fiiX3 + ¥3X1) = • 



(46) 
(47) 

(48) 
(49) 
(50) 



There are two independent solutions of this set of equations. One is given in eqs. (17) 

V2 



(19). A second solution is given by: 



A 




(Xi 



A 



B 



(51) 



where the factors A and B should obey to the relation 

hA 2 + f 2 B 2 + V3f 3 AB = . 



(52) 



At this level we assume that some additional input neglected so far in the analysis, such 
as for instance some choice of the soft supersymmetry breaking parameters, selects the 
first solution as the the lowest minimum of the scalar potential. 



It is important to note that the stability of the alignment in eqs. (17) and (18), under 



small perturbations can be proven. If one introduces small parameters in the VEV of the 
fields as follows 



{<Pi, 
A 




A 




A 



z 
1 



C, 



(53) 



it is only a matter of a simple algebra to show that, for small (z,Xi,x 2 ,yi,y 2 ), the only 
solution minimizing the scalar potential in the supersymmetric limit is indeed that given 
by (z,x 1 ,x 2 ,y 1 ,y 2 ) = (0,0,0,0,0). 
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Given the symmetry of the superpotential Wd, starting from the field configurations of 



eqs. (18)-p0[), (or ( 17 )-( 19 )) and acting on them with elements of the flavour symmetry 



group 5*4 x Z 4 , we can generate other minima of the scalar potential. Some of them are 

M ( n \ (<Pu) 




a : a 




(54) 



The new minima are however physically equivalent to those of the original set and it is 
not restrictive to analyze the model by choosing as local minimum the specific field con- 
figuration discussed so far. 

For the Froggatt-Nielsen field the non vanishing VEV is determined by the D-term 
associated with the U(1)fn symmetry. In fact, we regard this U(l) local symmetry. 

Since the field content displayed in Table 3 is anomalous, we assume the existence of 
additional heavy multiplets to cancel the anomaly. Here we do not need to specify these 
fields, but we must presume that their VEV's are all vanishing. The D-term in the potential 
is given by: 

V D = ^(M 2 FI -g FN \6\ 2 + ...) 2 (55) 

where g F N is the gauge coupling constant of U(1)fn and M FI is the Fayet-Iliopoulos term. 
The vanishing of Vd requires 

g FN \6\ 2 = M 2 FI . (56) 
Assuming that M FI / ' g F N is positive, this condition fixes the VEV of 9 and the value of t, 



given in eq. (21). Due to its particular origin the value of t is in principle unrelated to 



those of v and v' . 



6 Next-to-leading corrections 

We now study the set of subleading corrections to the superpotential that are essential 
to bring the model in agreement with the data. We start with the corrections to Wd that 
determines the vacuum alignment and then consider the modification to Wi and w u . The 
corrections can be classified according to an expansion in inverse powers of A and we will 
analyze the NLO terms to the LO results discussed before. 

6.1 NLO corrections to vacuum alignment 

The superpotential term Wd, linear in the driving fields ipf, x°, £° and <^°, is modified 
into: 

w d = w° d + Aw d . (57) 

where Awd is the NLO contribution, suppressed by one power of 1/A with respect to Wd- 
The corrective term Awd is given by the most general quartic, 5 4 x Z 4 x U(l) i?Ar-invariant 
polynomial linear in the driving fields, and can be obtained by inserting an additional 
flavon field in all the LO terms. The ^-charges prevent any addition of the flavons (pi 
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and xi & t NLO, while a factor of or <p v can be added to all the LO terms. The full 
expression of Awd is the following: 



\«=1 i=l i=l «=1 / 



(5? 



r £0 ^0 V)0 yO 1 

where £j, s, and Uj are coefficients and < J?", If", 7^ 1 , I i l > represent a basis of inde- 
pendent quartic invariants: 



n 



rv2 



r*2 



if 



L 

1 l 



if = &(<Pv(<Pv<P»)3) 

: (&i ( Pvh( ( PiXih) 



(59) 



(60) 



(61) 



v° 
J 3 



v° 



(62) 



((Xi^)s(^Xi)3) 

The new VEV configuration is obtained by imposing the vanishing of the first derivative 
of Wd + Awd with respect to the driving fields <^°, and Xi- We look for a solution 



that perturbs eqs. (17) and (18) to first order in the 1/A expansion: for all components 
of the flavons $ = (£„, (p„, <pi, xi)> we denote the shifted VEV's by 



($) = ($> LO + <5$ 



(63) 



where ($)lo are given by eqs. (17) and (18). 

After some straightforward algebra the results can be described as follows. In the 
neutrino sector the shifts 5£ v , 5<p v turn out to be proportional to the LO VEV's ($)lo 
and can be absorbed in a redefinition of the parameters C and D. Instead, in the charged 
lepton sector, the shifts 5<pi, 5\i have a non trivial structure, so that the LO texture is 
modified: 

Ml \ / $Xli 

{<*)= I A'A ( X i)= | (61) 

/ V B'A 



where A' and B' satisfy a relation similar to that in eq. (19) and the shifts and Sxii 
are proportional to v'vA, that are, in other words, suppressed by a factor v' with respect 
to the LO entries AA and BA, respectively. 
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6.2 NLO corrections to the mass matrices 

In this section we present the deviations from the LO mass matrices due to higher order 
operators and to the corrections to the VEV alignment. The NLO operators contributing 
to the lepton masses can be obtained by inserting in all possible ways or (p v in the LO 
operators and by extracting the S4 x Z4 x U(1)fn invariants. Insertions of one power of 
the flavons (pi or xi are forbidden by the invariance under the Z4 component of the flavour 
symmetry group. Now we list the terms that correct the mass matrices. 



6.2.1 Charged Leptons 



/ vt 2 vv't 2 vv't 2 \ 
mi = \ v't t Wd . 
\d' 1 / 



(66) 



The NLO operators contributing to the charged lepton masses originate from two 
sources: the first one is the introduction of the singlet in such a way that all the LO 
operators are multiplied by A and the resulting corrections are diagonal and of relative 
order v' with respect to the LO results; the second one consists in the insertion of the 
triplet (p u in the LO terms and the NLO superpotential is 



5wi 



A 3 A 2 



e c (l<fi(pi<p l 



v {2) e 2 

1 ^ e en 



yf) e 2 



A 2 A' 



A 2 A' 



(65) 



The charged lepton mass matrix is obtained by adding the contributions of this new set 



of operators, evaluated with the insertion of the LO VEV's of eqs. ( 17)-( 19 ), to those of 



the LO superpotential, eq. (14), evaluated with the NLO VEV's of eq. (64). By omitting 



all order one coefficients, we find 



In general each entry receives contributions from both the new operators in 5wi and from 
the LO terms of wi evaluated with NLO VEVs. Notice however that, at this order, the 23 
and 32 elements of mi are still vanishing. 

The unitary matrix that realizes the transformation to the physical basis where the 
product {m\mi) is diagonal at NLO is of the general form 

/ 1 V 12 v' V 13 v' \ 
Ui = -V 12 v' 1 , (67) 
V -V13V' 1 / 

where the coefficients are of 0(1). 
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6.2.2 Neutrinos 



The NLO operators contributing to the neutrino masses are given by 

+^( i /V c )(< /Wl/ ) + |((^) 2(v ^) 2 ) + (68) 

+j^((v c v c ) 3 ((p v (p v ) 3 ) + ^(i/ c i/V,)^ + h.c. . 

As we have already observed, the structure of the LO VEV's of the neutrino flavons is 
unchanged by the NLO corrections and the only possible modification to neutrino mass 
terms originates from the operators listed here. The Dirac mass matrix becomes 

yiD -y 2 C -y 2 C \ 

-y 2 C Vl D \v u . (69) 
-y 2 C Vl D / 

The Majorana mass matrix is also modified by the new operators: we note that the terms 
proportional to a± and to 04 can be absorbed by the LO term aD and we refer to all these 
contributions with the parameter a; the term proportional to 03 vanishes; the remaining 
term proportional to a 2 fills in new entries of Mjy and therefore it is the only genuine 
correction. The doublet component in the product of two ip v has a particular structure 
and its VEV preserves the subgroup generated by S: 

1 \ -,2*2 
















[! 







1 yv u + 




V 


1 


0/ 





((</w,) 2 > = {ja) C A ■ ( 7 °) 

Finally M N is modified into 

/ 2M + 2aD + 2a 2 C 2 -2bC -2bC \ 

M N = -2bC 3a 2 C 2 2M + 2aD - a 2 C 2 A . (71) 

V -2bC 2M + 2aD-a 2 C 2 3a 2 C 2 ) 

As we can easily see both Mjy and m u can be exactly diagonalized by the BM mixing, 
which represents the total contribution to lepton mixing of the neutrino sector, even after 
the inclusion of the NLO corrections. These corrections introduce also some terms in the 
mass eigenvalues of relative order v' with respect to the LO results, but they do not affect 
the type of the spectrum. 

Since the neutrino mass matrix is diagonalized by Ubm, the PMNS matrix can be 
written as 

U = ulU BM , (72) 

and therefore the corrections from Ui affect the neutrino mixing angles at NLO according 
to 

sin 2 12 = X - --L(V 12 + V 13 )v' 

sin 2 9 23 = X - (73) 
sin6» 13 = 7^(^12 - Vis)?/ • 
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ref. [21 

w L J 


ref. [3] 


sin 2 #12 


0.326+S-S5S [2<t1 

0.040 L J 


0.304+°-° 22 

— 0.016 


sin 2 #23 


0.45+2$ \2a] 

— u.uy l j 


0.50+g-gJ 

— u.uo 


sin 2 6>i3 


0.016 ±0.010 


n oin +0 - 016 

u.uiu_ 01 i 



Table 4: Results of two recent fits to the lepton mixing angles. 



By comparing these expressions with the current experimental values of the mixing 
angles in table |4j we see that, to correctly reproduce #12 we need a parameter v' of the 
order of the Cabibbo angle Ac- Moreover, barring cancellations of/among some the 
coefficients, also the reactor angle is corrected by a similar amount. 




0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 



Sin £12 



Figure 1: sin 2 #13 as a function of sin 2 6*12 is plotted, following eqs. (73 1. The plot is symmetric with 
respect sin 2 #12 = 0.5 and we report here only the left-part. The parameters Vij are treated as random 
complex numbers of absolute value between and 2, while \v'\ has been fixed at the indicative value of 
0.15. The gray bands represents the regions excluded by the experimental data [2]: the horizontal one 
corresponds to the 3<r-upper bound for sin 2 13 of 0.46 and the vertical ones to the region outside the 3a 
error range [0.26 — 0.37] for sin 2 Ow 



Any quantitative estimates are clearly affected by large uncertainties due to the pres- 
ence of unknown parameters of order one, as we can see in figure [TJ but in our model a 
value of #13 much smaller than the present upper bound would be unnatural. 

All this discussion works at the NLO, but we expect that at the NNLO the value of 
#23 will also be modified with deviations of about 0(\ 2 C ) at most. The next generation of 
experiments, in particular those exploiting a high intensity neutrino beam, will probably 
reduce the experimental error on #23 and the sensitivity on # 13 to few degrees. If no 
significant deviations from zero of #13 will be detected, our construction will be ruled out. 

A salient feature of our model is that, at NLO accuracy, the large corrections of 0(Xc) 
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only apply to 9 12 and #13 while #23 is unchanged at this order. As a correction of O(Xc) to 
#23 is hardly compatible with the present data (see Table 4) this feature is very crucial for 
the phenomenological success of our model. It is easy to see that this essential property 
depends on the selection in the neutrino sector of flavons £ v and <p v that transform as 1 
and 3 of S4, respectively. We have checked that if, for example, the singlet is replaced 
by a doublet ip v (and correspondingly the singlet driving field £° is replaced by a doublet 
ijj®), all other quantum numbers being the same, one can construct a variant of the model 
along similar lines, but in this case all the 3 mixing angles are corrected by terms of the 
same order. This confirms that a particular set of S4 breaking flavons is needed in order 
to preserve #23 from taking as large corrections as the other two mixing angles. 

7 Phenomenological implications 

We now develop a number of important physical consequences of our model and derive 
some predictions. We first consider the predicted spectrum and the effective mass m ee for 
neutrinoless double beta decay and then we discuss the constraints from leptogenesis. 



7.1 Light neutrino masses and neutrinoless double beta decay 

The light neutrino mass matrix, including the NLO corrections described in sections 
(4.2) and (6.2.2) is given by: 



(74) 



where m„, and m e J^ are given in eqs. (69), (71), (38), respectively. It is diagonalized 



by the BM unitary transformation and its complex eigenvalues are given by: 



m 2 



m 3 



2(F' + Y') - 
2(F' - Y') - 
-IF' + 



(y' + y 2 



2(F + Y + a 2 C 2 ) 

(2/ ~ Y 2 ) 2 
2{F-Y + a 2 C 2 ) 



K 
A 

v 2 

u 

X 



(75) 



,'2 



2(F - 2a 2 C 2 



K 
A ' 



where 



F' = M' + a'D 



Y' 



■V2b'C 



Y 



-V2y 2 C 



y' = y + yiD 



(76) 



We see that the LO expressions are recovered in the limit F' — Y' — y\ 2 — 0. By 



exploiting eqs. (76) we can study some observables like the effective 0z/2/5-decay mass, 
|m ee |, the lightest neutrino mass, mi, and the sum of the neutrino masses directly from 
the experimental data. We perform a numerical analysis, by treating all the LO, NLO 
and effective parameters as random complex numbers with absolute value between and 
3. In figure [2^ a), we plot | function of the lightest neutrino mass. The points 

displayed the figures correspond to the case of normal ordering of the neutrino masses, 
with a moderate hierarchy or a quasi degenerate spectrum. However, at variance with the 



18 



present Ov2/? bound: 



Majorana/GERDAIII 
CUORE 




10 -< 0.001 0.01 

m vl (eV) 

(a) \m ee \ vs. mi 



present 0v2y? bounds 




.. .GERDAIIT^ 

Majorana/GERDAIII 
t-CUORE 



Majorana phase a 2 \ 

(b) \m ee \ vs. a 2 i 



Figure 2: In figure (a), 



as a function of the lightest neutrino mass, mi, is plotted. The constraints 



which have been imposed in the plot are the experimental values at 3a for Am^ 



'so; 



and the mixing 



angles. All the parameters of the model are treated as random complex parameters. The present bound 
from the Heidelberg-Moscow experiment is shown in dark gray and the future sensitivity of CUORE (~ 15 
meV), Majorana and GERDA III (~ 20 meV), and GERDA II (~ 90 meV) experiments are represented 
by the horizontal dashed lines, while the future sensitivity of 0.2 eV of KATRIN experiment is shown by 
the vertical band. In figure (b), \m ee \ as a function of the physical Majorana phase olqx is shown. The 
experimental bounds are the same as in figure (a). 



results of the leading order, some solutions of our numerical simulation also reproduce an 
inverted hierarchy spectrum. The plot displays only the points corresponding to choices of 
the parameters reproducing Am^ (m , Am^ o/ and the mixing angles within a 3a interval. In 
dark gray, we plot the present bound from the Heidelberg-Moscow [18] experiment on \m ee \. 
The vertical band corresponds to the future sensitivity on the lightest neutrino mass of 
0.2 eV from the KATRIN experiment [19] and the horizontal lines to the future sensitivity 
of some 0z/2/5-decay experiments, that are 15 meV, 20 meV and 90 meV, respectively 
of CUORE [20], Majorana [21] /GERDA III [22] and GERDA II experiments. The figure 
suggests that a lower bound of about 10~ 3 eV holds for the lightest neutrino mass. Indeed, 
with the inclusion of the NLO corrections, from eq.(76) we see that m v \ can vanish if a 
cancellation between the NLO and LO contributions takes place. This however requires 
an additional fine tuning of the parameters which has been reproduced in our numerical 
analysis only partially and by very few points. Similarly the scatter plot indicates a lower 
bound for \m ee \ of about 0.1 meV. In figure [2^b), we plot | | clS db function of the 
physical Majorana phase. The maximum corresponds to a vanishing value of a.21 while 
the minimum to 021 = ±7r, as we expect from the general definition of |m ee |- In a non- 
negligible portion of the parameter space of the model, the predictions for |m ee | approach 
the future experimental sensitivity. Beyond the experiments already quoted we also recall 
SuperNEMO (50 meV) [33] and EXO (24 meV) [34]. 

In figure [3j we plot the sum of the neutrino masses as a function of the lightest neutrino 
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Figure 3: The sum of the neutrino masses as a function of the lightest neutrino mass mi. The vertical 
dashed line represents the future sensitivity of 0.2 eV from KATRIN experiment and the horizontal ones 
refer to the cosmological bounds (see the text for more details). 



mass, mi. The vertical band refers to the future sensitivity of KATRIN experiment, while 
the horizontal ones to the cosmological bounds [23] . There are typically five representative 
combinations of the cosmological data, which lead to increasingly stronger upper bounds 
on the sum of the neutrino masses: we are showing the two strongest ones. The first one 
at 0.60 eV corresponds to the combination of the Cosmic Microwave Background (CMB) 
anisotropy data (from WMAP 5y [24], Arcminute Cosmology Bolometer Array Receiver 
(ACBAR) [25], Very Small Array (VSA) [26], Cosmic Background Imager (CBI) [27] 
and BOOMERANG [28] experiments) plus the large-scale structure (LSS) information on 
galaxy clustering (from the Luminous Red Galaxies Sloan Digital Sky Survey (SDSS) [29]) 
plus the Hubble Space Telescope (HST) plus the luminosity distance SN-Ia data of [30] 
and finally plus the BAO data from [31]. The second line at 0.19 eV corresponds to all the 
previous data combined to the small scale primordial spectrum from Lyman-alpha (Lya) 
forest clouds [32]. This plot is typical for normal hierarchy or quasi degenerate spectrum. 
The only special feature is the lower bound on \mi\, which, as explained above, relies on 
a naturalness assumption. 



7.2 Leptogenesis 

In this section we briefly comment on the constraints on our model from the require- 
ment that it is possible to explain a sufficiently large baryon asymmetry through leptoge- 
nesis. The baryon asymmetry (see [35] for details) can be defined with respect the entropy 
of the universe as 

(8.75 ±0.23) x 10~ n , (77) 



_ n B - ng 

S 

where and s are the number densities of, respectively, baryons, antibaryons and 

entropy, a subscript implies "at present time", and the numerical value is from com- 
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bined microwave background and large scale structure data. When considering that the 
baryon asymmetry originates from leptogenesis, it is possible to link Yab and the lepton 
asymmetry parameter 6j as follows: 



'ab ^ C Y] 6i rji , (76 



i 

where C represents constants and r\i are the efficiency factors linked to the thermal equi- 
librium of the interactions of the RH neutrinos. We will not perform here a complete 
analysis of leptogenesis and in particular we will not evaluate the efficiency factors that, 
due to washout processes, reduce the final amount of lepton asymmetry. Rather, we 
formulate our constraint by requiring that the asymmetry parameters 6j are larger than 
about 1CT 7 (considering a quite usual value for the efficiency factor of about 0.1 and the 
supersymmetric context), by assuming that this value, after inclusion of the appropriate 
conversion factor, reproduces the observed baryon asymmetry. Since in our model RH 
neutrinos are typically very heavy, with masses above 10 12 GeV, we also assume that the 
correct framework to analyze leptogenesis is the unflavoured one. 

Out-of-equilibrium decays in the early universe of v c to lepton and Higgs doublets 
produce lepton asymmetries. In a basis in which the Majorana mass matrix for the RH 
neutrinos is diagonal and real, and where all basis-dependent matrices will be denoted by 
a hat, the lepton asymmetry parameters are 




where Y v = rh®/v u . Non- vanishing asymmetry parameters require that the off-diagonal 
entries of the product Y v Yj are complex and different from zero, while, at the LO, from 
eq. (23), we have 

t = U T BM Y U , Y v Yl = \y\ 2 l , (80) 

concluding that all the 6{ vanish at LO. 

At the NLO the asymmetry parameters still vanish. Indeed, as we have seen in the 
previous chapter, at NLO both Mn and Y v are still diagonalized by Ubm- A misalignment 
between and Y V YJ first occurs at the NNLO, when the VEV of the flavon field (p u 
acquires term of order v A /v' that do not preserve any more the symmetry generated by S. 
The first non-trivial contribution to the off-diagonal entries of Y v Yj is of order v 4 /v' and 
arises from the terms proportional to yi of the superpotential in eq. (68), when we take 
NNLO VEV's of the flavon tp v . In this case, considering that the lightest RH neutrino is 
z/f, leptogenesis in governed by e 3 : taking at this stage g ~ 1, we get 

1 v 8 

Considering the explicit expression of g we get a slightly different result. In this case the 
RH neutrinos are hierarchical and the function g scales as 



|m 3 | 2 ; " \M, 
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Figure 4: Contour lines of the asymmetry parameter, eq. (81 1 (continuous lines) and eq. (83) (dashed 
lines), in the plane (v,v'), for e 3 = (0.2,5) x 1CT 6 . 



In figure [4] we show the region of the parameter plane (v, v') compatible with the 
requirement e 3 ~ (0.2 5) x 10 -6 . Clearly, due to the scaling e 3 oc v 8 , the most sensitive 
parameter is v, which is bounded to a rather limited region and cannot be too small. 
Moreover, even a very large theoretical uncertainty, affecting our choice of the range e 3 ~ 
(0.2 -r 5) x 10 -6 , is considerably reduced when translated in terms of the parameter v. It is 
worth to stress that the region of v preferred by leptogenesis has a non- vanishing overlap 
with the values of v needed to explain the hierarchy in the charged lepton sector. 



8 Summary and conclusion 

We have constructed a see-saw model based on the flavour symmetry S4 x Z 4 x U(1)fn 
where the BM mixing is realized at the LO in a natural way. The hierarchy of charged 
lepton masses is obtained as a combined effect of the U(l) fn symmetry breaking measured 
by the parameter t and of the S4 x Z 4 breaking induced by v, proportional to the VEV's 
of ipi and xi- We have m^/m T — t ~ 0.06 and m e m T /m 2 ^ = v ~ 0.08. Within large 
uncertainties, a value v « (0.1 4- 0.2), not too far from the previous figure, is indicated by 
requiring that the observed baryon asymmetry in the Universe is explained by the amount 
of leptogenesis predicted by the model. 

Since exact BM mixing implies a value of tan # 12 which is excluded by the data, large 
corrections are needed. The dominant corrections to the BM mixing arise at NLO only 
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through the diagonalization of the charged lepton mass matrix. The shifts of the quantities 
sin 2 9 12 and sin #13 from the BM values are linear in the parameter v', proportional to the 
VEV's of ip u and £„, which is expected to be of the same order as v, but not necessarily 
too close, as v and v' are determined by two different sets of minimization equations. 
From the experimental value tan 2 12 = 0.45 ± 0.04, which is sizeably different than the 
BM value tan 2 12 = 1, we need v' ~ O(Xc) where ~ 0.23 is the Cabibbo angle. As 
in most models where the BM mixing is only corrected by the effect of charged lepton 
diagonalization, one also expects 9 13 ~ 0(\ c )- A value of 9 13 near the present bound 
would be a strong indication in favour of this mechanism and a hint that the closeness 
of the measured values of the mixing angles to the TB values may be purely an accident. 
In addition, a very important feature of our model is that the shift of sin 2 #23 from the 
maximal mixing value of 1/2 vanishes at NLO and is expected to be of O(Xq) at most. In 
our 64 model, this property is obtained by only allowing the breaking of £4 in the neutrino 
sector via flavons transforming as 1 and 3 (in particular with no doublets). 

In order to reproduce the experimental value of the small parameter r = Am 2 sun / Am 2 tm 
we need some amount of fine tuning. For instance, the RH neutrino Majorana mass M 
should be below the cutoff A (this is reminiscent of the fact that empirically M ~ M GUT 
rather than M ~ M P i anck ). The neutrino spectrum is mainly of the normal hierarchy type 
(or moderately degenerate), the smallest light neutrino mass and the 0vf3f3 parameter \m ee \ 
are expected to be larger than about 0.1 meV. The model is compatible with the observed 
amount of the baryon asymmetry in the Universe interpreted as an effect of leptogenesis. 
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Appendix A: the group S4 



We recall here the multiplication table for S4 and we list the Clebsch-Gordan coef- 
ficients in our basis. In the following we use ctj to indicate the elements of the first 
representation of the product and $ to indicate those of the second representation. 

We start with all the multiplication rules which include the 1-dimensional representa- 
tions: 

1 <S> Rep = Rep © 1 = Rep with Rep any representation 
1'® V = 1 ~ a (5 



V © 2 = 2 



a(3 2 
—afii 



V © 3 = 3' 





V © 3' = 3 



The multiplication rules with the 2-dimensional representation are the following ones: 

1 ~ (XlPi + Ct 2 P 2 



2©2 = l©l / ©2 



3 ©3' 



2 © 3' = 3 © 3' 



with 



with 



with 



1' ~ a±p 2 - 0£ 2 Pi 
2 ~ 



3' 



3' 



a 2 p 2 — 011P1 
ot\P 2 + ce 2 Pi 



ai(3 



\ 



^ot 2 p 3 - \aif3 2 
^a 2 p 2 - \a x p 3 j 



-a 2 Pi 
^a x p 3 + \a 2 p 2 
\ ^ot x p 2 + \a 2 p 3 ) 

( ~a 2 p! \ 

^aiAs + ^a 2 p 2 
\ ^«i/3 2 + \a 2 p 3 / 



Oil Pi 



\ 



^a 2 (3 3 - \a x p 2 
\ "§a 2 p 2 - \ ai p z J 
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The multiplication rules involving the 3- dimensional representations are: 

1 ~ aiPi + a 2 (3 3 + a 3 (3 2 

^{a 2 (3 2 + a 3 {3 3 ) 

a 3 (3 3 - a 2 (3 2 
ai(3 3 + a 3 Pi 
—ai/3 2 — a 2 (3i 

a 3 [3 2 - a 2 (3 3 
■V I a 2 Pi - ati(3 2 
a ± f3 3 - a 3 (3i 



3 © 3 = 3' © 3' = 1 © 2 © 3 © 3' with < 



( 1/ 



V ~ ai(3i + a 2 (3 3 + a 3 (3 2 



V3 



(a 2 (3 2 + a 3 (3 3 ) 



-aif3i + \ {a 2 f3 3 + a 3 /3 2 ) 



3 © 3' = V © 2 © 3 © 3' with < 



a 3 (3 2 - a 2 (3 3 
ai{3 3 - a 3 (3i 

a 3 p 3 - a 2 (3 2 
o I ati(3 3 + a 3 (3i 
—cti(3 2 — a 2 (3i 
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